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It is well known that non-central nuclear forces, such as the spin-orbital conpling and the ten¬ 
sor force, play important roles in understanding many interesting features of nuclear structures. 

However, their dynamical effects in nuclear reactions are poorly known since only the spin-averaged 
observables are normally studied both experimentally and theoretically. Realizing that spin-sensitive 
observables in nuclear reactions may carry useful information about the in-medium properties of 
non-central nuclear interactions, besides earlier studies using the time-dependent Hartree-Fock ap¬ 
proach to understand effects of spin-orbital coupling on the threshold energy and spin polarization 
in fusion reactions, some efforts have been made recently to explore dynamical effects of non-central 
nuclear forces in intermediate-energy heavy-ion collisions using transport models. The focuses of 
these studies have been on investigating signatures of the density and isospin dependence of the 
form factor in the spin-dependent single-nucleon potential. Interestingly, some useful probes were 
identified in the model studies while so far there is still no data to compare with. In this brief review, 
we summarize the main physics motivations as well as the recent progress in understanding the spin 
dynamics and identifying spin-sensitive observables in heavy-ion reactions at intermediate energies. 

We hope the interesting, important, and new physics potentials identified in the spin dynamics of 
heavy-ion collisions will stimulate more experimental work in this direction. 

PACS numbers: 25.70.-z, 21.30.Fe, 21.10.Hw 


I. INTRODUCTION 


Understanding novel features of the fundamental nu¬ 
clear forces and properties of strongly interacting matter 
under extreme conditions of density, temperature, spin, 
and isospin are among the main goals of nuclear physics. 
Heavy-ion collision (HIC) experiments play an important 
role in achieving these goals. Indeed, great achievements 
have been made using HICs at various beam energies 
from the sub-Colomb barrier to the highest energy avail¬ 
able at the large hadron collider. In particular, terrestrial 
experiments using intermediate-energy HICs have led to 
strong constraints on the equations of state of hadronic 
matter [1] and neutron-rich nucleonic matter mM- 


Theoretical studies have shown recently that some 
spin-sensitive observables of HICs can be used to ex¬ 
plore the in-medium properties of non-central nuclear 
forces. The spin-dependent nuclear interactions are im¬ 
portant for explaining several interesting features of nu¬ 
clear structure [1], such as the varying magic numbers 
and the shell evolution with the isospin asymmetry of fi¬ 
nite nuclei. However, the strength, density, and isospin 
dependence of the nuclear spin-orbit coupling are still 
uncertain (see Sec. HA). Moreover, the tensor force can 
modify the magic number of nuclei and is an impor¬ 
tant source of the nucleon-nucleon short-range correla¬ 
tion. The latter is related to many interesting phenom- 
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ena in nuclear physics (see Sec. HB). More studies on 
in-medium properties of the spin-orbit coupling and ten¬ 
sor force are thus vey much needed. HICs provide flexible 
ways of adjusting the conditions of the nuclear medium 
and may also lead to new spin-dependent phenomena. 
For example, the so-called “Spin Hall Effect ” ISHZ!, which 
affects the dynamics of spin-up and spin-down particles 
differently as a result of the spin-orbit coupling, is ex¬ 
pected to be a general feature in any spin transport pro¬ 
cess. It thus might be interesting to test if such phe¬ 
nomenon can also happen in HICs. 


Considerable efforts using the time-dependent 
Hartree-Fock (TDHF) model, the spin- and isospin- 
dependent Boltzmann-Uehling-Uhlenbeck (SIBUU) 
transport model, and the quantum molecular dynamics 
(QMD) model have been devoted to exploring the 
spin dynamics in HICs (see Sec. HI). Indeed, some 
interesting phenomena were found. For example, it was 
found that the inclusion of the spin-dependent nuclear 
interaction may affect the fusion threshold, generate 
the spin twist during the collision process, and lead 
to the spin splitting of nucleon collective flows (see 
Sec. HI). Future comparisons with relevant experimental 
data may help extract properties of the in-medium 
spin-dependent nuclear force. Here we review briefly the 
main physics motivations and recent findings of studying 
the spin-dependent dynamics and observables in low- 
and intermediate-energy HICs. A major goal of this 
article is to stimulate more experimental work in this 
direction. 
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II. SPIN-RELATED NUCLEAR FORCE 

Based on the one-boson-exchange picture [8], nuclear 
force can be understood by exchanging mesons between 
nucleons. Exchanging the scalar a meson and vector w 
meson leads to respectively the attractive and repulsive 
central nuclear force as well as the spin-orbit interac¬ 
tion, while exchanging the tt meson and p meson leads 
to respectively the long-range and short-range nuclear 
tensor force. Although in free space the bare nuclear 
force is well constrained by the nucleon-nucleon scatter¬ 
ing data, the in-medium nuclear interactions, especially 
the nuclear spin-orbit interaction and tensor force, are 
still quite uncertain. The in-medium nuclear interactions 
can be studied by using microscopic many-body theories 
or phenomenological models, such as the non-relativistic 
Skyrme-Hartree-Fock (SHF) model and the relativistic 
mean-field (RMF) model. In the following, we will dis¬ 
cuss the effective spin-dependent nuclear force based on 
the energy-density functional in the phenomenological 
approach. 


A. Nuclear spin-orbit interaction 

The nuclear spin-orbit interaction was first introduced 
to explain the magic numbers of nuclei lillD]. Nuclei 
with numbers of neutrons or protons equal to the magic 
numbers are more stable, and this reflects the special 
shell structure of a nucleus. Although even a simple har¬ 
monic potential leads to the shell structure of nucleon 
energy levels inside nuclei, the spin-orbit coupling is es¬ 
sential to reproduce the correct magic number. 

In the Skyrme interaction, the effective spin-orbit force 
between two nucleons at position fq and fq can be ex¬ 
pressed as m 

Vso = iWoiai + $ 2 ) -kx 6{fi - r 2 )k'. ( 1 ) 

In the above. Ho is the spin-orbit coupling constant, 
(Ti and S 2 are the pauli matrices for the two nucleons, 
k = —i(Vi — V 2)/2 is the relative momentum opera¬ 
tor acting on the right side with Vi and V 2 acting on 
the first and second nucleon, respectively, and k' is its 
complex conjugate acting on the left. From the conven¬ 
tional Hartree-Fock method, the spin-orbit single-particle 
potential can be obtained based on the above effective 
spin-orbit force 

K = ^9 ■ X (2) 

where 

W, = ^{Vp + Vp,) (3) 

is the form factor of the spin-orbit potential, with q = n 
or p being the isospin index and p being the nucleon 
number density. Taking the operator —iV as the mo¬ 
mentum p, the right-hand side of Eq. ([^ has the form 


of {f X p) ■ a with Wq playing the role of r, and this 
is why it is called the spin-orbit potential. By solving 
the Schrodinger equation with the single-nucleon Hamil¬ 
tonian 


with m* being the effective nucleon mass and Uq being 
the spin-independent potential, the single-nucleon spec¬ 
trum in a spherical closed-shell nucleus can be obtained. 

In the RMF model, Dirac equation is solved where the 
spin of nucleon is treated explicitly with nucleon wave 
functions for different spin states [l2] . The studies of SHF 
and RMF models on nuclear structure were reviewed in 
Ref. [T3], and here we compare the effective spin-orbit 
potentials from both the relativistic and non-relativistic 
approach. With non-relativistic expansion of the Dirac 
equation, the form factor of the nucleon effective spin- 
orbit potential in the RMF model can be expressed in 
the form of [HI [15] 


WrmF 


{2m 


1 


^eff'^P, 


( 5 ) 


where m is the nucleon mass and the coefficient Ceff is 
related to the coupling strength and mass of the scalar a 
meson and the vector w meson, i.e.. 


^ _ 9 I . al 

^eff — 9 9 ' 

ml ml 


( 6 ) 


The form factors of the spin-orbit potential in the SHF 
model (Eq. §) and the RMF model (Eq. (§) are dif¬ 
ferent. First of all, the spin-orbit coupling strength is a 
constant in the SHF model, but the effective coupling 
strength depends on the density in the RMF model. 
Implementing an additional density-dependent effective 
nucleon-nucleon spin-orbit interaction with a coupling 
constant Wi, the authors of Ref. m got addition con¬ 
tributions to the form factor as 

K = ^[cpy{p-p,) + {2 + c){2pqrVpq] 

+ ^Cp^~^{p- Pq)Vp, (7) 


with c mimicking the density dependence. The above 
form was tested in Ref. m in a semi-infinite nuclear 
matter with parameters fitted to the RMF interaction. 
It was found that the general features of the RMF model 
were then reproduced with this non-relativistic density- 
dependent spin-orbit interaction. Nevertheless, the den¬ 
sity dependence of the spin-orbit coupling is still largely 
unknown so far, and it is related to many interesting phe¬ 
nomena in nuclear structure studies [HHinj. Second, the 
spin-orbit couplings from the SHF and RMF approach 
have different isospin dependence, i.e., in the SHF ap¬ 
proach the spin-orbit coupling is stronger for nucleons 
of the same isospin, while in the RMF approach the cou¬ 
pling strength is the same for neutrons and protons. This 
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feature impacts descriptions of properties of neutron-rich 
nuclei, e.g., the kink in the evolution of the charge radii 
for lead isotopes. It was shown that the weak isospin 
dependence of the spin-orbit coupling in the RMF ap¬ 
proach can better explain the kink than the conventional 
SHF functional. However, if the form factor in the latter 
approach was modified to HOIISI] 

- Wn Wn 

w, = ^{i + x^)^Pg + ^^P,', {q^q') ( 8 ) 

a similar kink can be reproduced with Xw ~ 0.1 |21j . cor¬ 
responding to the case with very small Fock contribution 
of the spin-orbit interaction. Similar efforts were made 
by using a modified SHF functional to reproduce the 
isospin dependence of the spin-orbit field in semi-infinite 
nuclear matter with different neutron excesses |22] and 
in neutron-rich nuclei [ 23 j from a relativistic approach. 
In Ref. [ 23 ] , the isospin dependence of the spin-orbit cou¬ 
pling was compared in light drip line nuclei from the rela¬ 
tivistic mean field theory and the non-relativistic Skyrme 
model. Furthermore, it was observed that the commonly 
used Skyrme functional of the spin-orbit splitting over¬ 
estimated the central density and the spin-orbit splitting 
of neutron drops [ 25 ], calling for new functionals of the 
spin-orbit coupling. The proton energy splitting of /iii/2 
and ^7/2 outside the ^ = 50 closed shell increases with 
neutron excess, corresponding to the decreasing strength 
of the nuclear spin-orbit interaction | 26 |. The studies so 
far seem to favor a weak isospin dependence of the spin- 
orbit coupling. However, since the isospin dependence 
of the spin-orbit coupling, which is important in nuclear 
surfaces, is often coupled with its density dependence, it 
is still not well settled yet. 

Based on the above discussion, we proposed a general 
form of the form factor of the spin-orbit coupling by tak¬ 
ing both the density and isospin dependence into account 

= ^ + [q^q'). ( 9 ) 

The above form is artificially constructed and includes 
the main physics for simplicity purpose. In the above 
form, 7 is used to mimic the density dependence of the 
spin-orbit coupling while fixing its strength at saturation 
density po to be Wq. a and b are parameters to vary 
the isospin dependence of the spin-orbit coupling, with 
0 = 2 and 6=1 corresponding to the case of the stan¬ 
dard SHF approach and a = b corresponding to the case 
of the RMF approach. The values of 7, a, and 6 are 
still uncertain according to the above discussions. For 
the strength of the spin-orbit coupling Wq, efforts have 
been made to extract its information from ground-state 
properties of various nuclei. Recent studies have shown 
that the spin-orbit coupling and the tensor force, which 
will be discussed in the next subsection, should be con¬ 
sidered simultaneously to describe the spin-orbit split¬ 
ting and single-nucleon spectra of nuclei. Based on the 
Skyrme functional and taking the uncertainties of the 


tensor force into account, the strength of the spin-orbit 
coupling is approximately 80—150 MeV fm®, from fitting 
the properties of light to heavy nuclei . 

The single-nucleon Hamiltonian of Eq. 0 is adequate 
to describe the ground-state properties of spherical closed 
shell nuclei. For open shell nuclei, one needs to consider 
an additional spin-dependent potential using the spin- 
current density J from Eq. 0 

= •(/+/,). (10) 


J is actually the vector component of the spin-current 
density tensor see, e.g.. Ref. m for more detailed 
discussions. Eor deformed nuclei, not only the time-even 
potentials (Eqs. 0 and @) but also the time-odd po¬ 
tentials should be considered m- Starting from the ef¬ 
fective nucleon-nucleon spin-orbit interaction (Eq. 0 ) 
and taking Eq. 0 into consideration, the general form 
of the time-even and time-odd spin-dependent potentials 
can be written as 


U~ = [V • (aJ, + 6 /,,)] 

+ ^ {aVpq + bVpg,)-{pxa), ( 11 ) 

^ P- ^ + ^V)] 

- ^ x{ajq + bjq,)],{q^q') 

( 12 ) 


where p = —iV is the momentum operator, and s and 
j are spin density and current density, respectively. We 
note that the time-odd potentials play an important role 
in the dynamics of heavy-ion reactions, which will be 
discussed in Sec. nD 


B. Nuclear tensor force 

The first strong evidence of the nuclear tensor force 
is from studying properties of deuterons. In the non- 
relativistic approach, the nuclear tensor force between 
two nucleons at position fq and Iq is often expressed with 
the tensor operator written as 

o 0(^1 (' 5^2 -rO 7- Mo^ 

012 — 3 - ^ -(cti • CT2)) ( 13 ) 

where r = fq — fq is the relative position vector. One 
thus sees that whether the tensor force is attractive or 
repulsive depends on the relative direction between the 
spin and the relative position vector, i.e., 512 > 0 for 
(fi(2) parallel to r and 5i2 < 0 for ffi(2) perpendicular to 
r. 

The tensor term in the effective nuclear interaction was 
first included in the Skyrme force | 31 | . but afterwards it 
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was neglected due to its complex form. Recently it at¬ 
tracted renewed interests. It has been found by Otsuka 
et al. that the nuclear tensor force may affect the shell 
structure or even modify the magic number of nuclei [51- 
151] . The combined effects of the spin-orbit coupling and 
the nuclear tensor force sometimes hamper our under¬ 
standing on both of them [3^ [35] . Based on the random 
phase approximation, effects of the tensor force on the 
multipole response of magic nuclei have been studied, 
and a large effect on the magnetic dipole states was ob¬ 
served m- Besides, it was proposed that the spin-isospin 
excitation of finite nuclei may serve as a useful observ¬ 
able to assess the strength of the tensor force [3H1133] . 
The existence of the tensor force may also open a shell 
gap for large neutron numbers, having a consequent im¬ 
plication for the synthesis of neutron-rich superheavy el¬ 
ements [40) . 

Although the nuclear tensor force has no effect on 
the equation of state of spin-saturated nuclear matter 
based on the studies at mean-field level, it affects the 
properties of nuclear matter from many-body calcula¬ 
tion beyond the mean field. It was found that the re¬ 
pulsive central force and the tensor force are two im¬ 
portant sources of nucleon-nucleon short-range correla¬ 
tion |41] . The high-momentum tail of nucleon distri¬ 
bution in nuclear matter as well as in finite nuclei was 
observed even at zero temperature based on these stud¬ 
ies |33H35] . Great efforts have been made in measuring 
the short-range nucleon-nucleon correlation and extract¬ 
ing the ratio of nucleons in the high-momentum tail ex¬ 
perimentally [351(50] and theoretically [511 [53] , see, e.g.. 
Ref. |53j for a review. In particular, it was found that 
the neutron-proton correlation is much stronger than the 
correlation between neutron-neutron and proton-proton 
pairs [35050] , and this is mainly due to the nuclear tensor 
force. The isospin dependence of the short-range corre¬ 
lation can lead to interesting consequences, such as the 
reduction of the kinetic contribution to the nuclear sym¬ 
metry energy |54H56j compared to the free Fermi gas sce¬ 
nario. Since the symmetry energy at saturation density 
is constrained to be around 30 MeV from many analyses, 
see, e.g.. Ref. m, the isospin-dependent short-range cor¬ 
relation effectively increases the potential contribution 
to the symmetry energy and thus the symmetry poten¬ 
tial effect, which may lead to enhanced isospin effects in 
intermediate-energy HICs [55] [53] . 


III. SPIN IN NUCLEAR REACTIONS 

The spin hall effect was first predicted by M.I. 
Dyakonov and V.I. Perel in 1971 [ 313 , while the term 
’’Spin Hall Effect” was named by Hirsch in 1999 [5]. Con¬ 
sidering the transport of spin-up and spin-down particles 
with spin-orbit coupling, i.e., = —L ■ a with L being 

the angular momentum and a being the particle spin, the 
spin-up (spin-down) particles tend to turn left (right) to 
couple with the angular momentum and lower the energy. 



Spin Hall effect 



FIG. 1: Cartoon illustrating the spin hall effect. Taken from 

Ref. [5]. 



leading to the splitting of final observables for different 
spin states, as shown in Fig.[2 In this section, we will dis¬ 
cuss similar effects in low- and intermediate-energy HICs 
based on the framework of TDHF, BUU, and QMD, with 
more complicated forms of spin-orbit coupling from spin- 
dependent nuclear interactions. 


A. TDHF model study 

The mean-field dynamics of nucleons in the TDHF 
model is described by 


i 


dt 




h(f>i , 


(14) 


where (j)i is the wave function of the ith nucleon and the 
single-nucleon Hamiltonian is given by 


h(j)i = 



(15) 


with E being the energy functional of the nuclear sys¬ 
tem from Hartree-Fock calculation. The single-nucleon 
Hamiltonian is generally a function of nucleon number 
density p, spin density s, current density j, spin-current 
density J, and so on, and their definitions in terms of the 
nucleon wave function are 


P = 

i 

S = ^^(?i*(cr|d|cr')0i, (17) 

i ( 7 ,( 7 ' 

7 = ( 18 ) 

i 

J = X (cr|CT|cr'), (19) 


with (o’ld’lcr') being the pauli matrix element. Numer¬ 
ically, these densities can be calculated on the coor¬ 
dinate space grid and Eq. can be solved with a 

fixed time step. Eor more details, we refer the reader 
to Refs. [601 El]- The TDHE framework works well for 
low-energy heavy-ion reactions and in studing resonances 
dynamics. 
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FIG. 4: Spin excitation in central reactions from 

TDHF study with both time-even and time-odd terms. Taken 
from Ref. [5S| . 


FIG. 2: The relation between the out-going energy and in¬ 
coming energy in O^® -I- O^® reaction from TDHF study. The 
filled (open) circles are results with (without) spin-orbit force, 
and the solid, dashed, and dotted curves represent results 
from three different Skyrme forces. Taken from Ref. m- 


In the old calculations, the single-nucleon Hamiltonian 
was generally calculated from the SHF model without 
the spin-orbit interaction and time-odd terms. The spin- 
orbit force was first introduced to the TDHF framework 
in Refs. [5^64) . It is interesting to see that the spin-orbit 
force enhances the dissipation in the fusion reaction and 
transforms the relative motions of the two nuclei into the 
internal excitations. The fusion threshold energy in 0^®-|- 
reaction is increased by about a factor of 2 [5^ IS5] . 
as shown in Fig. with three different parameterizations 
of Skyrme force. The fusion cross section obtained from 
TDHF calculation was increased after including the spin- 
orbit force [5^ . 



FIG. 3: Center-of-mass energy evolution in central O^® -I- O^® 
reaction with (solid lines) and without (dashed lines) time- 
odd contributions from TDHF calculation. The time window 
of the reaction process is from about 50 fm/c to 120 fm/c. 
Taken from Ref. |65| . 


With only time-even contribution of the spin-orbit in¬ 
teraction, i.e., the spin-orbit potential (Eq. (§) and the 


potential with spin-current density J (Eq. (10)), spuri¬ 
ous spin twist can be generated in a free moving nucleus, 
as a result of spin-orbit coupling. Obviously, this phe¬ 
nomena is not reasonable as it depends on the reference 
frame. Considering that all kinds of collision geometry 
can be realized in HICs, the time-odd terms were further 
introduced in the TDHF calculation in Refs. [5SJ [SS] to 
satisfy the invariance under Galilei transformations. It 
is seen from Fig. that there is no such spurious spin 
and the kinetic energy is a constant before 50 fm/c when 
the nuclei are moving freely, as a result of suppression 
effect on the time-even terms from the time-odd terms. 
During the reaction process, the real spin twist appears 
due to the overwhelming effect of the time-odd terms on 
the time-even terms, as shown in Fig. At the end 
of the reaction, the energy of outgoing nuclei is smaller 
with the time-odd terms as shown in Fig. indicating a 
stronger dissipation. Besides the spin excitation, it was 
found that the fusion description was further improved 
with the time-odd terms and the spin-current pseudoten¬ 
sor contribution |66] . A more detailed study on this topic 
was done recently m, where it was found that the dissi¬ 
pation is dominated by the time-even contribution of the 
spin-orbit force at lower energies but by the time-odd 
terms at higher energies. 


(a) Protons 



(b) Neutrons 



-SLy5+T 
■■ SLy4d+Stancu 
-SkM*+Slancu 
-SV-tls 


Mass number 


Mass number 


FIG. 5: (Color online) Ratios of the spin mean field from the 
spin-current density representing the tensor force to that from 
the spin-orbit force for protons (left) and neutrons (right) as 
a function of the mass of reaction nuclei. Different parame¬ 
terizations of the Skyrme force are used in the study. Taken 
from Ref. |55| . 
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Besides the spin-orbit force, the additional contribu¬ 
tion of the spin-current density J was introduced in the 
TDHF calculation representing the contribution from the 
tensor force in Ref. jSS] • It was found that the dissipation 
effect from the tensor force is small compared with that 
from the spin-orbit force |68[l69j . However, the spin mean 
field can be enhanced with the nuclear tensor force, and 
the enhancement becomes important with the increasing 
mass of the colliding nuclei, as shown in Fig. depend¬ 
ing on the parameterization of the Skyrme force. A more 
complete study by including the full Skyrme functional 
as well as the tensor force in the TDHF calculation was 
done very recently in Ref. [70]. It was found that the 
Skyrme tensor force has non-negligible effects on low- 
energy heavy-ion dynamics and the fusion threshold en¬ 
ergy. 


B. BUU model study 


spin degree of freedom are described by 


dr p Wo f p\^ ^ ^ . 

Jt = m + —(,«) -x(aVp, + i.Vft.) 

^ V X (asg + 


dp 

dt 

da 

dt 


Po 


s—odd 


-VUg - - VUg 


Wo ) [{aV Pg -I- bVpq') X i^x d 

kFo (^ ) [V X {ajg + bjg>)] X a. 


( 20 ) 

( 21 ) 


( 22 ) 


One sees that the three degrees of freedom couple with 
each other. The number density p, the spin density s, 
the current density j, and the spin-current density J are 
calculated from test particle method inii7ii7i|. Since 
the mixing of the long-range Fock contribution and the 
spin interaction is a complex problem, the momentum 
dependence is not included in the spin-independent po¬ 
tential Uq for the moment. In addition, the spin of nu¬ 
cleons are randomized after nucleon-nucleon scatterings, 
by approximately taking the spin flip effect into consid¬ 
eration IZilTS]. 


The TDHF model works well in low-energy HICs, while 
the particle emission and nucleon-nucleon scattering are 
still lacking. To describe these effects in intermediate- 
energy HICs, BUU models and QMD models are suit¬ 
able candidates. In the BUU framework, the Boltzmann 
equation is solved with test particle method [iDlig. In 
the previous studies, an isospin-dependent BUU (IBUU) 
transport model has been used to describe the isospin 
dynamics in intermediate-energy HICs (2. Recently, the 
spin degree of freedom of nucleons and the spin-orbit in¬ 
teraction were incorporated in the IBUU model, and the 
new model is dubbed as the spin- and isospin-dependent 
BUU (SIBUU) model jTSHZZj- In this section, we summa¬ 
rize the main results published originally in Refs. [73ll77] . 

In the SIBUU model, each nucleon is assigned ran¬ 
domly a unit vector representing the expectation value 
of its spin. In this way, the spin projection of each nu¬ 
cleon at arbitrary direction can be easily calculated. In 
the transport simulation, z direction is set as the beam 
momentum and x direction is for the impact parameter. 
Since the total angular momentum in non-central HICs 
is in the y direction perpendicular to the reaction plane, 
i.e., X — o — z plane, it is reasonable to study the spin 
polarization in y direction. We thus determine the nu¬ 
cleons with spin projection on +y {—y) direction as the 
spin-up (spin-down) nucleons. 


Considering the general form of the time-even and 
time-odd spin-dependent potentials in Eqs. (III and (121, 
the time evolutions of the coordinate, momentum, and 



FIG. 6: (Color online) Time evolution of contours of reduced 
number density p/po, V component of the spin density Sy, 
X component of the number density gradient (Vp)^;, and y 
component of the curl of the current density (V x j)y in non¬ 
central Au-|-Au collision at the beam energy of 50 MeV. Taken 
from Ref. [73]. 


The time evolution of relevant density contours from 
SIBUU calculation are displayed in Fig.[^ The gradient 
of number density Vp and the curl of the current den¬ 
sity j show the strength of the time-even and time-odd 
spin-dependent potential, respectively, and both of them 
are closely related to the evolution of the number density 
shown in the first row of Fig. [^ The nucleon spin tends 
to be oarallel to p x Vp from the time-even potential 
(Eq. (Ill), while it tends to be parallel to V x j from the 
time-odd potential (Eq. (iHl))- The contributions from 
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the time-even and time-odd potentials are opposite to 
each other. One sees that before the two nuclei touch 
each other there is no spin polarization as a result of the 
cancellation of the time-even and time-odd potentials, 
consistent with the foundings from TDHF studies. Dur¬ 
ing the collision process, the participant is polarized in 
the -\-y direction, i.e., in the direction of the total an¬ 
gular momentum, following the preference direction of 
the time-odd potential. It is seen that the direction of 
the spin polarization is consistent with that in Fig.|^from 
TDHF calculation with both time-even and time-odd po¬ 
tentials. 


Au+Au@50MeV, b = 8 fm 



CD 

< 

c? 


FIG. 7: (Color online) Transverse flow of spin-up and spin- 
down nucleons (left) and spin up-down differential transverse 
flow (right) with different strength of the spin-orbit coupling 
in non-central Au-|-Au collisions at the beam energy of 50 
MeV. Taken from Ref. [73] . 


Transverse flow is one of the most important ob¬ 
servables for extracting the equation of state of pro¬ 
duced matter and studying the nuclear interaction in 
HICs m HU). The left panel of Fig. displays the 
transverse flow of spin-up and spin-down nucleons as a 
function of reduced rapidity ■ We note that the 

target (projectile) nucleus is in the +x (—x) direction in 
Fig.§ which is different from the conventional initializa¬ 
tion, leading to the negative slope of the transverse flow. 
However, this doesn’t prevent the reader from seeing the 
obvious splitting of transverse flow between spin-up and 
spin-down nucleons. With a detailed orientation analy¬ 
sis, one can find that again the time-odd potential domi¬ 
nates the effect, giving the spin-up (spin-down) nucleons 
an attractive (repulsive) potential. This can be under¬ 
stood in a naive picture that the spin-up (spin-down) 
nucleons parallel (antiparallel) to the direction of total 
angular momentum and thus feel an attractive (repul¬ 
sive) potential. One can further define the spin up-down 


differential transverse flow as follows 


Fud{yr') 


1 

^iVr) 


NiVr) 

i=l 


(23) 


where ai is 1 for spin-up nucleons and —1 for spin-down 
nucleons, and N(yr) is the number of nucleons at rapidity 
yr- The above spin up-down differential transverse flow 
largely cancels the effect from the spin-independent nu¬ 
clear interaction while preserves the information of the 
spin-dependent potential. Indeed, the slope of Fud in¬ 
creases with increasing spin-orbit coupling constant, in¬ 
dicating that it is a good probe of nuclear spin-dependent 
interaction. 



FIG. 8: (Golor online) Dependence of the spin np-down dif¬ 
ferential transverse flow on the beam energy (left) and the 
impact parameter (right) in non-central Au-|-Au collisions. 
Taken from Ref. [7i| . 

The spin up-down differential transverse flow was fur¬ 
ther analyzed in detail in Ref. [Ti]. Figure [^displays the 
dependence of F^d on the beam energy and the centrality. 
At higher beam energies, the angular momentum is larger 
while the nucleon-nucleon scattering is more violate, with 
the former enhancing the spin-dependent potential while 
the latter washing out part of the information of spin dy¬ 
namics. The competition leads to a maximum slope of 
Fud at the beam energy of about 100 MeV, as shown in 
the left panel of Fig. Since the spin-dependent poten¬ 
tial is related to the density gradient and is thus a surface 
effect, the slope of Fud increases with the increasing value 
of the impact parameter, as shown in the right panel of 

Fig.m 

In the neutron-rich collision system where the relevant 
neutron densities are larger than proton densities, the dif¬ 
ference of the spin up-down differential transverse flow of 
neutrons and protons can be useful to probe the isospin- 
dependence of the spin-orbit coupling in HICs. The anal¬ 
ysis was carried out with a stronger isospin-like coupling 
(a = 2, 6 = 1) and a stronger isospin-unlike coupling 













FIG. 9: (Color online) Spin up-down differential transverse 
flow of neutrons and protons at different beam energies and 
with two typical isospin dependence of the spin-orbit coupling. 
Taken from Ref. [7T| . 



FIG. 10: (Color online) Slope parameter of the spin up-down 
differential transverse flow F' without (left) and with (right) 
high transverse momentum cut from different density depen¬ 
dence of the spin-orbit coupling at different beam energies. 
Taken from Ref. [Til¬ 


ing spin-orbit coupling strength, compared to the case 
with a constant one. In this way, the strength and the 
density dependence of the spin-orbit coupling can be dis¬ 
entangled. 


(a = 1, 5 = 2), and the resulting F^d were calculated 
at different beam energies shown in Fig. A stronger 
isospin-like spin-orbit coupling, which is exactly the case 
of SHF interaction, leads to a larger F^d for neutrons 
than for protons, while a stronger isospin-unlike coupling 
gives opposite predictions or similar Fud for neutrons and 
protons. The effect is appreciable from beam energy 50 
MeV to 200 MeV, while the beam energy of 100 MeV is 
the optimized one due to the largest magnitude of F^d- 

The density dependence of the spin-orbit coupling has 
bothered many nuclear physicists and hampered the un¬ 
derstanding of nuclear spin-orbit interaction in nuclear 
structure studies. Since HICs have the advantage of con¬ 
structing the system with designed density, isospin, and 
momentum current, it might be helpful in extracting use¬ 
ful information of the density dependence of the spin- 
orbit coupling. As is known, nucleons of high transverse 
momentum (pt) emit early from the high-density phase 
in HICs, and the density of the high-density phase in¬ 
creases with increasing beam energy. This feature can 
be used to extract the density dependence of the spin- 
orbit coupling, as illustrated in Fig. Without high¬ 
er cut, the slope of F^d can hardly be distinguished as 
shown in the left panel of Fig. because nucleon emis¬ 
sion from low-density phase, which is similar at different 
beam energies, dominates the results. With high-pT cut, 
the slope of F^d is smaller at lower collision energies but 
larger at higher collision energies from a linearly increas- 



FIG. 11: (Color online) Transverse momentum dependence of 
the elliptic flow of mid-rapidity nucleons in peripheral Au-|-Au 
collisions at different beam energies. Taken from Ref. ESI- 

In non-central HICs, the azimuthal distribution of 
emitted nucleons can always be expressed as 

d^N d^N 

-b 2v2{yr,PT)cos{2(t)) + ..] (24) 

with (j) = [jpy/Px) being the azimuthal angle, and 

vx = (cos((/>)) and t >2 = (cos(20)) are called the directed 
flow and elliptic flow, respectively. The directed flow 
is similar to the transverse flow but it depends on the 
flow angle rather than magnitude. The elliptic flow is 
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positive at lower energies, negative at intermediate ener¬ 
gies, and becomes positive again at higher energies. The 
positive elliptic flow means more particles move in-plane 
than out-of-plane as a results of hydrodynamics, while 
the negative elliptic flow is a result of the squeeze-out ef¬ 
fect on the expansion of participant matter by the spec¬ 
tator nucleons [T]. Despite of the complicated dynamics, 
the elliptic flow serves as a useful probe of the proper¬ 
ties of nuclear matter formed in HICs and the nuclear 
interaction. The transverse momentum dependence of 
V 2 of spin-up and spin-down nucleons at mid-rapidity is 
displayed in Fig. El Except for the different behaviors 
of V 2 at different beam energies, the large elliptic flow 
of spin-up nucleons than spin-down nucleons is observed, 
especially at higher transverse momentum as a result of 
the stronger spin-orbit coupling for energetic nucleons. 
At the energy range considered, a more attractive mean- 
field potential leads to a larger V 2 in peripheral HICs, 
consistent with the effect of spin-dependent potential on 
the spin splitting of transverse flow discussed above. 

The above observables are for free nucleons. Exper¬ 
imentally it is easier to detect charged particles rather 
than neutrons, leading to difficulties of measuring the 
spin splitting of transverse flows for protons and neutrons 
and identifying the isospin dependence of the spin-orbit 
coupling. Of course the spin measurement is another 
challenge which will be discussed in the next section. 
Once the corresponding detectors are set up, the spin 
splitting of observables for charged light clusters may be 
more easily measured. For transport models with point¬ 
like particles, the dynamical coalescence approach has 
been shown to be successful in studying the hadroniza- 
tion in relativistic HICs [STJ [55] and light cluster forma¬ 
tion in intermediate-energy HICs jSSl [13]. In this ap¬ 
proach, the probability for nucleons to form a light clus¬ 
ter is proportional to the nucleon Wigner function of the 
light cluster [531 HI] ; and the proportional constant is the 
statistical factor determined by the spin-isospin degener¬ 
acy. For example, with explicitly knowing the isospin of 
nucleons, the statistical factor for a neutron and a pro¬ 
ton to form a deuteron is 3/8, while that for one neutron 
and two protons to form a ^He is 1/12. Since now the 
spin of each nucleon is also explicitly known, the dynam¬ 
ical coalescence can be further improved by considering 
the antisymmetrization of the product of spin and isospin 
wave function. For example, the statistical factor for a 
spin-up neutron and a spin-up proton to form a spin-up 
deuteron is 1/2, while for a spin-up neutron, a spin-up 
proton, and a spin-down proton to form a spin-up ^He is 
1/2. This improvement has been applied to study spin 
splitting observables for deuterons, tritons, and ^He [55] . 
It has been checked that after spin average the results 
reproduce those without explicit spin treatment. 

Figure displays the spin splitting of the directed 
flows for deuterons, tritons, and ^He in non-central 
Au-|-Au collisions at the beam energy of 100 MeV. The 
directed flow of spin-down clusters is larger than that of 
spin-up ones. The spin splitting of the directed flow is 



FIG. 12: (Color online) Directed flow of deuteron, triton, and 
^He of different spin states in non-central Au-|-Au collisions 
at the beam energy of 100 MeV. represents the spin state 
projecting on the y direction perpendicular to the reaction 
plane. Taken from Ref. [85] . 

largest for deuterons due to its large spin quantum num¬ 
ber, i.e.. S' = 1. The spin splitting observables of tritons 
and ^He might be more easily measurable for extracting 
the isospin dependence of the spin-orbit coupling exper¬ 
imentally. 



FIG. 13: (Color online) Elliptic flow of of spin-up and spin- 
down deuterons in non-central Au-|-Au collisions at the beam 
energy of 100 MeV. Sz represents the spin state projecting on 
the y direction perpendicular to the reaction plane. Taken 
from Ref. [85] . 

The elliptic flow of different spin states of deuterons 
has been illustrated in Fig. [T^ in non-central Au-|-Au 
collisions at the beam energy of 100 MeV. It is seen that 
the elliptic flow of spin-down deuterons is more negative 
at mid-rapidity, but is slightly positive at large rapidity, 
indicating an obvious spin splitting even taking the sta¬ 
tistical error into account. Again, the magnitude of the 
V 2 as well as its spin splitting for deuterons is larger than 
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that of free nucleons according to Ref. m, and might 
serve as a better spin-dependent observable. 



FIG. 14: (Color online) Spin up-down differential transverse 
flow from full Skyrme calculation using MSLO force without 
(left) and with tensor force (right). Taken from Ref. | 85 |. 


Further preliminary calculation with full Skyrme func¬ 
tional has been done. A standard Skyrme functional 
with MSLO parameterization [SS] has been used in the 
calculation. The detailed derivation and expression of 
full Skyrme functional with both time-even and time-odd 
terms can be found in Refs. ETIIMIISO]- The resulting 
spin up-down differential transverse flow is shown in the 
left panel of Fig. 14 One can see the similar sensitiv¬ 
ity of Fud to the spin-orbit coupling strength although 
the magnitude is a little smaller, compared to the result 
shown in Fig. where only the spin-orbit coupling is ap¬ 
plied. To investigate the effect of nuclear tensor force on 
the spin dynamics of intermediate-energy HICs, a zero- 
range tensor force of the form 


vt{r) = ■ k'){a2 ■ k') - ■ a2)k'‘^]6{r) 

+ ^(0[3(cti ■ k){a2 ■ k) - (di • CT2)fc^]} 

-I- to[3((Ji • k')5{r){a2 ■ k) - (cti • a 2 )k' ■ 5{r)k] 

(25) 


is incorporated to the full Skyrme transport model calcu¬ 
lation, where r = ri — r 2 is the relative coordinate, k and 
k' are the relative momentum operator and its complex 
conjugate, respectively, and te and to are the triplet-even 
and triplet-odd strength parameter. The energy density 
function derived from the above tensor force can be found 
in Refs. Eziiin], where the corresponding terms (such as 
the spin-current density J) are non-negligible only when 
local spin polarization is produced. The resulting spin 
up-down differential transverse flow is shown in the right 
panel of Fig. It is seen that the slope of F^d is not 


very sensitive to the values of tg or to unless extremely 
large coupling constant is used. This feature is consistent 
with TDHF study where the spin dynamics is dominated 
by the spin-orbit coupling. However, one would expect 
that with spin-polarized beam or target, the tensor force 
effect can be much enhanced. 


C. QMD model study 


In the QMD framework, the Wigner function of each 
nucleon is treated as Gaussian wave packet in both co¬ 
ordinate and momentum space m 1^ . and the two- 
nucleon interaction is related to the effective two-body 
interaction and the overlap of their wave functions. The 
equation of motion in the QMD model is given by the 
semiclassical canonical equation, i.e., 


df 

dt 

dp 

dt 


VpiL, 

-VrH, 


(26) 


where f and p are respectively the central coordinate and 
momentum of the wave packet, and FI is the Hamiltonian 
of the system including the kinetic and potential energy. 

In a recent study, the nuclear spin-orbit interaction was 
incorporated to the ultra-relativistic QMD (UrQMD) 
model. The potential energy contribution of the spin- 
orbit interaction is expressed as [SS] 


Us 



(27) 


where the spin-dependent potential Us consists of the 
time-even and time-odd contribution written as 


ur^ = -^(pv-j+p„v-4), (28) 
^odd ^ _:^[,-..(vxj) 

+ Sn ■ (V X jn) -|- Sp • (V X jp)], (29) 


where Wq represents the spin-orbit coupling strength, 
and p, s, j, and J are the number, spin, current, and 
spin-current densities, which can be calculated from lo¬ 
cal Wigner function of the nucleon [M]. 

The spin dynamics was analyzed based on the above 
framework. Similar spin splittings of the directed flow 
and the elliptic flow were observed in non-central Au-|-Au 
collisions at the beam energy of 150 MeV, as shown in 
Fig. lA It was argued that the net spin-dependent po¬ 
tential is attractive for spin-up protons and repulsive for 
spin-down protons, leading to a larger directed flow for 
spin-down protons than spin-up protons. The spin split¬ 
ting of pT-integrated elliptic flow was found to be small 
and only visible in peripheral collisions, and it was found 
that V 2 for spin-down protons is slightly larger than that 
for spin-up ones. Since the conventional initial direction 
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Au+Au, E|^|^=150 MeV/nucleon, Free protons 



FIG. 15: (Color online) The directed flow (upper panels) and 
elliptic flow (lower panels) for spin-up and spin-down protons 
in non-central Au-FAu collisions at the beam energy of 150 
MeV from QMD calculation. Taken from Ref. [89]. 


of the target and projectile is used as shown in Fig. 3 of 
Ref. [53] , the spin-up (spin-down) nucleons correspond to 
the spin-down (spin-up) ones in the SIBUU study [751 - 
[77]. Although the spin splitting of final collective flow is 
a robust phenomenon in both models, further studies are 
needed to understand the relative sign of the splitting. 

Analysis was further done on the beam energy depen¬ 
dence of the flow splitting. Defining K„p and Kdown as the 
slope parameter of the directed flow of spin-up and spin- 
down protons, the slope difference is shown to increase 
with increasing impact parameter, as shown in the up¬ 
per panel of Fig.[^ qualitatively consistent with SIBUU 
studies. In non-central Au-|-Au collisions, it was found 
that the slope difference first increases then decreases 
with increasing beam energy, and the maximum differ¬ 
ence appears at the beam energy of 150 MeV, similar 
to the finding in the SIBUU model where the maximum 
slope of the spin up-down differential transverse flow ap¬ 
pears at the beam energy of about 100 MeV. 

It was further emphasized in Ref. [53] that the spin 
averaged flow results do not change after including the 
spin-orbit interaction, as a result of the cancellation of 
spin-up and spin-down nucleons. In addition, the spin 
splitting of the flow slope caused by the spin-orbit in¬ 
teraction is comparable to the isospin splitting caused 
by the nuclear symmetry energy, especially for neutrons. 
These findings are all consistent with the observations in 

SIBUU studies [73HZ2]- 


IV. EXPERIMENTAL STATUS 

Due to the difficulties of spin measurement in HIC 
experiments, the main focus in the past is mainly on 
the spin-averaged observables, so that the information 


Au+Au, Free protons 




FIG. 16: (Golor online) The difference of the slope of directed 
flow between spin-up and spin-down protons as a function 
of impact parameter (upper panel) and beam energy (lower 
panel). The dashed line in the upper panel is a linear fit, 
while the inset in the lower panel shows the relative difference. 
Taken from Ref. [89]. 


of spin dynamics is neglected. Thanks to the great ef¬ 
forts made by experimental nuclear physicists, the mea¬ 
surements of the spin of free nucleons and light clus¬ 
ters now become possible. Although the detailed ex¬ 
perimental status will be presented in another topic re¬ 
view of this issue, here we’d like to briefly mention 
two related experiments that might be relevant in an¬ 
alyzing the probes discussed above. One of them is 
the spin-polarized beam which can be produced through 
pick-up or removal reactions at Rikagaku Kenkyusho 
(RIKEN) [33] [3T], Gesellschaft fiir Schwerionenforschung 
mbH (GSI) [33], the National Superconducting Gyclotron 
Laboratory (NSGL) [93], and the Grand Accelerateur Na¬ 
tional d’lons Lourds (GANIL) [35H33] . It is expected 
that the effects of spin dynamics with spin-polarized 
beam will be much enhanced, providing a better system 
for extracting the information of the spin-dependent nu¬ 
clear force, especially the nuclear tensor force. For the 
spin-excitation state of heavy clusters, the spin polar¬ 
ization and alignment can be measured via the angu¬ 
lar distribution of its 7 or /3 decay, see, e.g.. Ref. m 
for a review. Making use of the analyzing power of a 
nucleus might be the most promising way of identify- 
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ing the spin of free nucleons or light clusters experimen¬ 
tally. The analyzing power indicates the left-right scat¬ 
tering asymmetry of an incident polarized nucleon on 
the target nucleus. The spin-dependent scattering is a 
result of the interference of electromagnetic interaction 
and hadronic force [98l, and the spin flip is observed 
between not only charged-charged scatterings but also 
charged-neutral scatterings. It is noteworthy that at cer¬ 
tain energies and scattering angles the analyzing power 
can be as large as 100% [99]. Experimental efforts are 
thus encouraged by using the selected nucleus as a ’de¬ 
tector’ whose analyzing power is known in prior. In this 
way the spin of corresponding particles can be measured 
and the probes discussed in the previous sections can be 
analyzed. 

V. SUMMARY 

In summary, we outlined the major physics motiva¬ 
tions of investigating the in-medium spin-dependent nu¬ 
clear interactions, i.e., the spin-orbit interaction and the 
nuclear tensor force, and summarized some recent efforts 
in exploring the spin dynamics in low- and intermediate- 
energy heavy-ion collisions. In particular, the studies 
on the strength, the density, and the isospin dependence 
of the spin-orbit interaction as well as the short-range 
correlation induced by the tensor force are highlighted. 
In TDHF studies, it has been found that the spin-orbit 
interaction can enhance the dissipation in low-energy 
heavy-ion reactions and increase the fusion threshold. 
Incorporating both the time-even and time-odd contri¬ 
bution of the spin-orbit interaction can lead to nontrivial 
spin polarization, while the tensor force slightly enhances 


the spin field compared to the spin-orbit interaction. In 
the studies using HICs at intermediate energies, both 
the spin- and isospin-dependent BUU model and QMD 
model predict the spin splitting of the nucleon collective 
flow, which seems to be a robust phenomenon. In the 
BUU model studies, efforts have been made in extract¬ 
ing the isospin dependence of the spin-orbit coupling and 
disentangle its strength and density dependence. Prelim¬ 
inary results on spin splitting of observables related to 
the light clusters and those from full Skyrme calculation 
with nuclear tensor force have also been discussed in the 
BUU model studies. We hope the findings summarized in 
this review will soon stimulate more experimental work 
in this direction. 
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